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FOREWORD 


Part of this report was written at the Department of Mathematics, 
Monash University, Clayton, Victoria 3168, Australia, and part during 
tenure of an NRC-NASA research associateship at the Jet Propulsion 
Laboratory. 

A shorter version of this report will be published in Physics 
Letters. 


ABSTRACT 


Theoretical modelling of the contraction of the primitive proto- 
saturnian cloud, using ideas of supersonic turbulent convection, suggests 
that each of Saturn'*/ inner moons, excepting Rhea, condensed above the ice- 
point of water and consists primarily of hydrous magnesium silicates. It 
is predicted that Voyager 1 may find that the satellite mean densities 
steadily increase towards the planet and that the rocky moons are irregular 
in shape. 
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1 . INTRODUCTION 


The recent spate of sightings [1] of suspected new satellites of 
Saturn observed during the rare ring-plane crossing of this planet in 
March 1980 has aroused a fresh interest in the origin of this planet and 
its inner family of moons and rings. This interest is particularly 
intense in view of the forthcoming fly-by of Voyager 1, scheduled 
November 12, 1980. Very little is known about the inner moons. Photo- 
metric evidence suggests that they have frost covered surfaces and low 
mean densities typical of a mostly icy composition [2]. In this report 
I present the results of detailed calculations of a possible mechanism 
for the formation of the moons and rings, which is based on ideas of 
supersonic turbulent convection applied to the original Laplacian 
hypothesis [3]. The calculations, which are a sequel to the work of 
Prentice and ter Haar [4], suggest that all of the moons interior to the 
orbit of Rhea, as well as the particles of the F ring discovered by 
Pioneer 11 [5] and the E ring [6], consist mostly of hydrated silicates. 
The A, B, C rings are mostly icy. These results can be reconciled with 
the existing observational data if the moons are irregular in shape. 

The orbital radii of the inner family of moons (Mimas through 
Rhea) form a roughly geometric sequence of the form Rp/R^^+l = 1 .30 t 0.07, 

where Rq denotes the orbit of Rhea and n = 0,1,2 The moon 1979 S 1 

discovered by Pioneer 11 [5], as well as the Fountain-Larson object 1966 
S 2 [7], lie amongst a suspected faint stream of satellites whose mean 
orbital radius of 2.52 R^ (Rg = Saturn's equatorial radius 60,000 km) 
is close to the next term in the geometric sequence after Mimas. The 
regular spacings of these orbits, including their near circularity and 
coplanarity and common sense of motion around Saturn, are features which 


1 


are also seen in the regular satellite systems of Jupiter and Uranus 

and the Sun's own planetary system. It is natural therefore to assume 

that some common mechanism was responsible for the formation of all systems [8], 

2. THE MODERN LAPLACIAN THEORY 

One possible mechanism which seems fairly accurately to account 
fnr the mass and chemical composition of the Galilean moons and the 
planetary system is a variant of the original nebula hypothesis due to 
Laplace (9], In this modern Laplacian theory [3, 4) the satellites and 
planets condensed from a concentric system of orbiting gaseous rings that 
were shed by the rotating turbulent gaseous envelope which gravitationally 
contracted to form the central parent body. 

In the case of Saturn, this gaseous envelooe was captured by a rock 
+ ice planetary core of mass 'v 10 (Mg = mass of Earth) which had accreted 

Q 

on the mean circular orbit of a protosolar gas ring of mass ^ 10 Mg and 
temperature 60 K. This gas ring had, in turn, been shed by the turbulent 
protosun at the orbit of Saturn during the Sun's own primordial contraction. 

Since the temperature of this gas ring is much less than the H 2 O ice point 
(150 K), one would expect the gas to be very heavily depleted of H 2 O, 
including admixed rocky condensate, A large fraction of the uncondensed 
gases, consisting mostly of H and He woulo thermally escape from the solar 
system prior to their capture by the core when the orotosun passed through 
its overluminous phase [3,X0]. The model of the primitive p^'otosaturnian 
cloud would thus consist of a central dense planetary core of mass 10 Mg 
surrounded by a tenuous rotating gaseous envelope of mass 85-90 Mg. The 
heavy element abundance of the envelope is expected to be much less than 
the solar value Zg = 0.018 [li] and its spin axis would be parallel to that 
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of the orbital motion. The present mass of Saturn is Mg » 95.2 Mg. 

It is interesting to note that Anderson et al. [12] have established 
an entirely similar picture from completely different considerations. 

These authors have computed the present structure of Saturn using 
gravitational moment data returned by Pioneer 11. They have concluded 
that H 2 O and possibly NHg and CH^ are primarily confined to the vicinity 
of a 15 - 20 core, if they are present in solar proportions to rock. 

Consider now the evolution of a tenuous protosaturnian envelope 
whose size initially exceeds the orbital radius of Rhea. Pollack et al. 

[ 13 ] showed that such a cloud would be fully convective, but unable to exist 
in a state of hydrostatic equilibrium since its internal energy requirement 
exceeds the available gravitational energy. Instead the cloud would collapse 
hydrodynamically over a period of a few days to about the orbit of Mimas, 
allowing no time to dispose of its excess spin angular momentum efficiently 
through the orderly process of mass extrusion and shedding [ 14 ]. To over- 
come this problem, which also manifests itself in the study of both the 
protosolar and protojovian envelopes, Prentice [ 3 ] has suggested that a 
supersonic turbulent stress arising from the motions of overshooting con- 
vective elements may help stabilize the envelope against free collapse. 

The turbulent convective stress is buoyancy driven and in a non- 
rotating cloud has the form - gpGM(r)/r where P is the gas density, 

M(r) the mass interior to radius r and 6 is the gravitation constant [ 15 ]. 

The turbulence parameter g'v 0.1 is assumed to be a constant throughout the 
cloud. The turbulent stress equals some 10 times the normal gas pressure 
pRT/p (R = gas constant, T = temperature, y = molecular weight) in the outer 


convective layers. It drastically lowers the cloud's moment-of-inertia 
coefficient f [16], thereby increasing the absolute store of gravitational 
potential energy. The stress also leads to the development of a very 
dense shell of non-turbulent gas above the photosurface. When rotation 
is included, this outer shell evolves into a belt-like structure of mass 
m at the equator, which is discontinuously abandoned by the contracting 
cloud when the centrifugal force overcomes the gravitational force. The 
net result is that the rotating cloud of mass M sheds a concentric system 
of orbiting gaseous rings whose mean orbital radii satisfy the equation 

Wi ■ ” * 

If the contraction is uniform, meaning that both m/M and f are constant, 
then the sequence is geometric. Eq. (1) can be inverted to estimate 
empirically the masses of the gaseous rings and hence compute the expected 
masses of the condensing system of moons, given the value of f and the 
envelope composition. In the case of the Gal 'lean moon system, a best fit 
to the observed masses occurs with the choice f = 0,02, for a gas of solar 
composition [4]. 

3. CALCULATIONS 

Applying eq. (1) to the protosaturnian cloud, with M = Mg, f = 0.02 

27 

and Rp/Rp+i = 1*3, we obtain a mean ring mass m = 1.6 x 10 g. If this gas 

24 

is of solar composition the mass of rock condensate per ring is 7 x 10 g 

25 

whilst the rock + HgO ice component is 1.9 x 10 g, using Cameron's 
abundances [11]. We see from Table 1, however, that both these masses are 
much larger than the observed values of any of the regular moons. If Rhea 
is mostly icy, the shortfall in expected mass corresponds to a factor of 9. 
This suggests that the heavy element abundance Zg of the saturnian envelope 
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was depleted relative to the solar value Zg ty the same factor, yielding 

Zg ■ 0.002. 

Now applying eq. (1) specifically M Dfone's orbit with Zg « 0.002 
a'.id Rq/Rj * 1.40 we obtain a rock mass of 1.0 x 10 g. This precisely 
coincides with the observed mass of this moon suggesting that Dione is, in 
fact, rocky, The steady decline in the masses of the moons interior to 
Dione can also be accommodated in terms of a rocky composition if much of 
the rock remains suspended in the gas ring as a fine dust unable to 
settle onto the mean orbit before the ring dispersed, taking with it the 
finer material [4]. The presence of the gas ring is essential both for 
focusing attains onto the mean circular Keplerian orbit of the gas ring 
[3,17] as well as for the efficient accretion of the settled material by 
the growing sateltite embryo [18]. 

At Dione’s orbit the time scale t^^^ for the settling of grains of 
radius a » 0.01 cm and density pg s 3 g cm is 150 yr whilst the dispersal 
time scale, based on a simple thermal evaporation argument [4], is 

Here T^^ = 275 K is the gaseous ring temperature (see table 1), Moving in- 
wards from Dione, che rings are warmer (T^^ » so the grains may be 
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dustier and smaller with a correspondingly longer settling time t^g^ = 1/t . 

At the orbit of Mimas t^g^ = 2600 yr if a = 10"^ cm, while tg^^p = 30 yr, 
Hence, only a very small fraction of the rocky condensate may have settled 
out onto the mean orbit by the time the gas had dispersed. Thermal 
stirring of the dense inner gas rings, whose density is shown in table 1, 
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would also destroy the orbital angular velocity structure in the ring, 
causing the gas and Its Imprisoned dust to spread away from the mean 
accretlonal oroit R^. We suggest, therefore, that Felbelman’s E ring 
formed from the slowly settling dust which failed to become Incorporated 
into the rocky moons [19]. 

The temperatures of the gaseous rings at their time of detachment 
t^ from the contracting globe can be computed by evolving a turbulent 
pc-lytropic model of mass M ■* Mg through the dimensions of the regular 
satellite system C203, The results of these computations for a cloud 
having H mass fraction X « 0.75 and ^ 0.002 are shown in fig, 1. This 

V 

figure shows the variation of the black-body surface temperature T^ of the 
cloud with equatorial size R^. During uniform contraction, T^ * ^"^^^ph^^1 
where Tp^j(R^.) Is the photosurface temperature at the Initial cloud radius 
R^. « 10 Rg. The values T^ * orbital positions of the satellites 

and mid-point positions of the A, B, F rings, Including Cassini's division, 
are shown In table 1, along with the expected composition of the condensate 
that would form, derived from the equilibrium condensation sequence data of 
Lewis [21], 

4. PREDICTED CHEMISTRY, DENSITIES AND SHAPES OF THE SATELLITES 

The simplified polytropic calculations indicate that Rhea condensed 
below the water vapour ice point. This moon should therefore consist of 3655 
by weight rock and 64% ice, with a mean density close to 1.3 g cm” [22J. 

The gravitational energy released during the accretion of Rhea's outer mantle 
would be sufficient to melt the icy mix if this moon's mass exceeds 2.5x10^^ g. 
It follow' that the surface of Rhea may be nearly spherical. Closer to Saturn 
the moon?: are predicted to be increasingly rockier with a mean density of 
order 2.4 g cm'^, similar to that of asteroidal rock [21], Peale et al. [23], 
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using tidal considerations, have also predicted that Mimas and probably 
Olone may be rocky. Thermodynamic data of hydrous silicates [24] suggest 
that the masses of the moons are just sufficient to liberate some of the 
water of hydration from the outer mantles during accretion. This steam 
'..ould form a torus around the circula*' orbit of the satellite which would 
subsequently recondense on Its surface to give It a frosty appearance, 
when the temperatures everywhere began to drop. In this connection 1979 
SI and Mimas should be darker and denser than the other moons since the 
mineral phase tremollte Is less hydrous than serpentine. 

The masses of the rocky moons are far too small for them to have melted 
during accretion and assumed a spherical form. We cannot, therefore, over- 
look the possibility that they may be Irregular in shape. There is an 
apparent discrepancy between the radiometric estimate of Dione's diameter 
[23 and the lunar occultation value [25], taken when the satellite was at a 
different orbital phase, namely 68® before eastern elongation [26]. Dione 
also has a non-s1nusoidal light curve [27]. These data suggest that this 
moon is a rocky elongated body, having a semi-major axis of about 575 km 
directed towards Saturn and minor axes of 425 km perpendicular to this. 
Tethys has a larger visual magnitude than Dione, despite its smaller mass. 
Combining this fact with Elliot's uniform disk measurements [25] and Tethys' 
uniform lignt curve, we deduce that this moon. If rocky, is a bright oblate 
spheroid having semi-major axes of about 450 km in the orbital plane and a 
minor axis of 300 km normal to this. Enceladus has a much lower photo- 
metric density than Mimas. It also has an unusual light curve [28], This 
moon may therefore be elongated like Dione, but of higher visual albedo and 
having a longest semi-axis of 300 km. A schematic Illustration of the pre- 
dicted satellite shapes and sizes is shown in Fig. 2. 
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5. FORMATION OF SATURN'S RINGS 


During the slow final stages of the cloud's contraction, the degree 
of turbulence measured by the parameter e is progressively reduced, 
commencing at radius ' 2.67 Rg, in a manner such that the surface 
temperature Tg follows the locus shown in fig. 1. = 300 K is the 

final chosen photosurface temperature at radius R = 1.37 R^., corresponding 
to a fully rotating cloud having polar radius equal to the present value 
0.912 Rg [5], As B declines, the mass of the shed rings and their spacing 
R^ "%+l rapidly decreases. With R* = 2.67 R^ the cloud sheds one final 
ring at orbital radius 2.52 Rg, equal to that of 1 979 S 1, then remains 
rotationally stable until the radius 2.42 Rg is reached. Bflow this point 
the cloud can remain stable only if f starts to increase. The outer non- 
turbulent layers, however, cannot be braked during contraction and will be 
stripped away creating a gaseous Keplerian disc in the equatorial plane. 

In our model it is the material which condenses from this disc which forms 
Saturn's rings [29]. 


Initially the condensate in the disc is rock. After a time t 


HgO’ 


typically about 3x10 yr, the temperature at each point falls below the 
H^O ice-point and ice begins to precipitate out of the gas. Owing to the 
thermal evaporation, however, the fraction of gas ())|^ g, shown in table 1, 
which survives to time t^ g decreases rapidly with distance from Saturn. 

The calculations thus show that whilst the B and C rings may be mostly icy, 
the fraction of ice to rock rapidly decreases through the region of the A 
ring. The particles of the F ring and other suspected nearby moonlets should 
therefore all be rocky [30]. Finally, according to this theory Titan and 
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the outer moons of Saturn are captured bodies which originally condensed at 
the same orbit as Saturn but failed to be accreted Into Its central core. 

As such, they should consist of 32% by weight rock, 57% H 2 O Ice and 11% NH 3 . 

6 . CONCLUSIONS AND ACKNOWLEDGEMENTS 

Voyager 1 will test the predictions of this theory In November, 1980. 
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Fig. 1. Surface temperature of the contracting protosaturnian cloud of mass 
Ms = 5.69 X 10^^ g plotted against equatorial size Rg, which is 
measured in units of the present value Rs = 6 x lO^ cm. The cloud 
is assumed to be in a state of supersonic turbulent convective 
equilibrium and uniform internal rotation, with centrifugal force 
balancing the gravitational force at the equator. The broken lines 
show the vapour condensation temperatures of H2O ice and the hydrous 
mineral serpentine Mgs [Si20s] (0H)4 at the cloud's equator. These 
lines were computed using the equilibrium condensation data of 
Lewis [22] applied to a gas of heavily depleted metal abundance 
Z$ = 0.002. The positions of Saturn's rings are shown at the foot 
of the diagram. 
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Fig. 2. Schematic illustration of pre-Voyager predictions of the shapes, 
sizes, and densities of Saturn's inner satellites; sizes refer to 
semi -axes (km). 
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